Short-Distance Correlation Properties of the Lieb-Liniger System and Momentum 
Distributions of Trapped One-Dimensional Atomic Gases 
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We derive exact closed-form expressions for the first few terms of the short-distance Taylor expansion 
of the one-body correlation function of the Lieb-Liniger gas. As an intermediate result we obtain 
the high-p asymptotics of the momentum distribution of both free and harmonically trapped atoms 
and show that it obeys a universal law for all values of the interaction strength. We discuss the 
ways to observe the predicted momentum distributions experimentally, regarding them as a sensitive 
identifier for the Tonks- Girardeau regime of strong correlations. 



Introduction. Even though the correlation functions 
for the Lieb-Liniger gas of ^-interacting one-dimensional 
bosons [1] have been an object of intense research in 
the Integrable Systems community since the late 70s [2], 
the full closed-form expressions are known only in the 
Tonks-Girardeau limit of infinitely strong interactions 
[3]. While the scaling properties of the long-range [4] 
asymptotics of the correlation functions can be derived 
from Haldane's theory of quantum liquids [5], Confor- 
mal Field Theory [6], and Quantum Inverse Scattering 
method [2,7], virtually nothing is known about short- 
range one-body correlations at finite coupling strength 
[8] . One of the goals of this paper is to extend the exist- 
ing knowledge in this direction. 

It is known that while for weak interactions the Lieb- 
Liniger system is well-described by the mean-field the- 
ory, the opposite, Tonks-Girardeau regime of infinitely 
strong interactions [9,10] constitutes a strongly correlated 
system dual to a free Fermi gas. In experiments with 
one-dimensional atomic gases [11,12] the one-body mo- 
mentum distribution of the gas, along with the density 
profiles [13] and phase fluctuations [14-16], can readily 
help to distinguish between the two quantum regimes. 
In the Tonks-Girardeau limit, the momentum distribu- 
tion for both free and harmonically confined gases was 
investigated by several authors [3,17-19]. In this paper, 
we address the question of the momentum distribution 
in the intermediate, in between mean-field and Tonks- 
Girardeau, regime, as more realistic from the experimen- 
tal point of view. 

System of interest. Consider a one-dimcnsional gas of 
N ^-interacting bosons confined in a length L box with 
periodic boundary conditions. The Hamiltonian of the 
system reads 
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where m is the atomic mass, and gm is the one- 
dimensional coupling constant, whose expression for real 
atomic traps is given in [17]. This Hamiltonian can be 
diagonalized via Bethe ansatz [1]. At zero temperature, 
the energy of the system is given through 



EN = n^eM 

' 2m ^" 
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where the dimensionless parameter 7 — 2/n]aiD| is in- 
versely proportional to the one-dimensional gas param- 
eter n [fliDl, ri is the one-dimensional number density of 
particles, aiD = —2h^/mgii) is the one-dimensional scat- 
tering length introduced in [17], and the function 6(7) is 
given by the solution of Lieb-Liniger system of equations 
[1] : it is tabulated in [20] . Note the asymptotic behavior 
of 6(7) (first computed in [1]): 
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where 7 — s- corresponds to the mean-field or Thomas- 
Fermi regime, whereas 7 — s- cx) corresponds to the Tonks- 
Girardeau regime. 

High-p momentum distribution. Our first object of in- 
terest is the high-p asymptotics of the one-body momen- 
tum distribution in the ground state. To evaluate it, we 
need two mathematical facts, (a) and (b): 

(a) The presence of the delta-function interactions in 
the Hamiltonian (1) implies that its eighenfunctions un- 
dergo, at the point of contact of any two particles i and 
J, a kink in the derivative proportional to the value of 
the eigenfunction at this point [21]: 
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^{zi, . . . , Zi, . . . , Zj, . . . , zn) 
= ^'(zi, . . . , Zji, . . . , Zji, Zn) (5) 
X {1 - \zj,\/am+e(\zj,\; {Zji})} 

where Zji = {zi + Zj)/2 and zji = zj — Zi 
arc the center-of-mass and relative coordinates of 
the ij pair of particles, respectively, and {2ji} = 
{Zji, z\, . . . , Zi-\, Zi+i, . . . , Zj-i, Zj+i, . . . , Zn} de- 
notes a set consisting of the center-of-mass coordinate 
of the i-th and j-ih. particles and the coordinates of all 
the other particles. 

(b) Imagine that a periodic function f{z), defined on 
the interval [—L/2, +L/2\, has a singularity of the form 
/(z) = \z — zq\"-F{z), where F{z) is a regular function, 
Q: > —1 and a ^ 0, 2, 4 . . .. Then the leading term in 
the asymptotics of the Fourier transform of / reads [22] 
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where k = {2n /L)s and s is an integer. For multiple 
singular points of the same order, the full asymptotics is 
the sum of the corresponding partial asymptotics of the 
form (6). 

Let us evaluate, using (5) and (6), the momentum 

representation of the ground state wavefunction of the 
Hamiltonian (1) with respect to the first particle: 
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Here pi = {2nU/L) s, where s is an integer. 

Let us now turn to the one-body momentum distri- 
bution per se. After a lengthy but straightforward [23] 
calculation it takes the form 

I-+L/2 r+L/2 

w{p) = / dz2... dzNl'i'ip, Z2, ■ ■ ■ , zn)\'^ 

J-L/2 J-L/2 
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where p^izx, Z2, z[, Zg) is the two-body density matrix, 
normalized as J^l^l2 dzi J^l^l2 dz2 P2{zi, Z2; ^1, ^2) = 1, 
and w{p) is the momentum distribution, normalized as 

E^r-ooM27rWi) = l. 

The expression (8) involves the two-body density ma- 
trix whose form is unknown for a finite system. However, 
an elegant thermodynamic limit formula for p2(0, 0, 0, 0) 
does exist due to Gangardt and Shlyapnikov [15], who 
derived it using the Hellmann-Feynman theorem [24]: 
^2(0, 0, 0, 0) = e'{j). We are now ready to give a 
closed-form thermodynamic limit expression for the high- 
p asymptotics of the one-body momentum distribution 
for one-dimensional J-interacting bosons in a box with 
periodic boundary conditions: 
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where W{p) = {L/2TTfi) w{p) is normalized as 
J-^ dpW{p) = 1. Notice that this asymptotics is uni- 
versally described by a law for all values of the cou- 
pling strength 7. (Note that for 7 00, this law was 
predicted in [19].) Formula (9) is the first of the two 
principal results of our paper. 

Harmonically trapped ID gas: momentum distribu- 
tion. To evaluate the high-p asymptotics of the momen- 
tum distribution of atoms confined in a harmonic oscil- 
lator potential, we are going to employ the local den- 
sity approximation (LDA). It is based on an intuitive, 
but hard to justify rigorously assumption that in the 
thermodynamic limit the short-range correlation prop- 
erties of a trapped gas are indistinct from the ones of 

z — z' 

a uniform gas of the same local density: a{z, z') « 
cr„(z — z' I n((z + z')/2), where cr(z, z') and (t„(z — z' \ n) 
are one-body density matrix of the trapped gas and one- 
body density matrix of a density n uniform gas respec- 
tively, both normalized to the respective number of par- 
ticles, n{z) is the density profile of the trapped gas, I is 
the typical length on which the density changes. (Tested 
against the exact results on correlation function of the 
trapped Tonks-Giradeau gas [25] this assumption can be 
shown to lead to an exact prediction for the value of 
the coefficient in front of the {z — z')'^ term in the Tay- 
lor expansion around z' = —z k, 0.) From this ansatz 
it immediately follows that the high-p asymptotics of the 
momentum distribution (sensitive to the short-range cor- 
relations only) is given by the spatial average of the uni- 
form case expression (9) over the density profile of the 
atomic cloud. The density profiles themselves can also 
be obtained using LDA ( [13], Eqns. 21 and 22, where 
the governing parameter 77 should be replaced by 2/7^p, 
see below), and this is the method we used. The final re- 
sult is presented in Fig.l. There the dimensionless coeffi- 
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cient n = lim|j,|^oo W {p)p / {hn")^ in front of the high-p 
(8) asymptotics of the momentum distribution is plotted as 
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a function of the interaction strength parameter 7*^ in 
the center of the cloud; 7° in turn depends on of the 
experimental parameters, such as the number of parti- 
cles TV, the coupling constant giD, and the longitudinal 
trap frequency u, through a system of implicit equations 
[13]. Here nP is the density in the center of the trap. 
To establish a link to the experimental parameters we 
also present a plot for the Thomas-Fermi (weak interac- 
tions) prediction for 7 in the center of the atomic cloud, 
7^F = (8/32/3)(7Vma2cj/?i)-2/3, as a fimction of 7°. 

In the limiting, Thomas-Fermi and Tonks-Girardeau 
regimes the momentum distribution is given by 
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FIG. 1. Dimcnsionlcss coefficient Q,Ho{'~f) in front of the 
leading term of the high-p asymptotics of the momentum dis- 
tribution of harmonically trapped atoms, as a function of 
the interaction strength 7° in the center of the cloud. Also 
shown is the directly related to the experimental parameters 
Thomas-Fermi estimate 7?.p = (8/3^''''*)(iVma?Qt^'/fi)-2/3 

the interaction strength in the center, as a function of 7*^. 
Both momentum and momentum distribution are measured 

in units related to the density in the center n", that can be 
expressed through 7" using 7" = 2/n"aiD. 

Short-range expansion for the correlation Junction. Let 
us now redirect our attention to the ground state one- 
body correlation function 

g,{z) = {¥{z)i-{Q)), (11) 

and in particular to its Taylor expansion around zero: 



gi{z)/n = 1 + ^Cj' 
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In the limit of infinitely strong interactions 7 
expansion is known to all orders [3]: 
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Our goal now is to obtain the first few (through the order 
|zp) coefficients of the expansion (12) for an arbitrary 
interaction strength 7. 

The knowledge of the momentum distribution (9) is 
crucial for determining the Ci and C3 coefficients. Let 
us look at the relation between the momentum distri- 
bution and the correlation function, where the former 
is simply the Fourier transform of the latter: W{p) = 
(27rnn)-i j^^dze-^P^/^'-giiz). Since the leading term 
in the asymptotics of W(p) is 1/p'^ we may conclude, 
using the Fourier analysis theorem (6), that the lowest 
odd power in the short-range expansion of the correla- 
tion function gi{z) is \z\^, and therefore the \z\ term is 
absent from the expansion: 



ci =0 



(14) 



Furthermore, the theorem (6) allows one to deduce the 
coefficient C3 from the momentum distribution (9): 



C3 = ^ 7^ e'(7) . 



(15) 



To obtain the coefficient C2, we employ the Hellmann- 
Feynman theorem [24] again. Let a Hamiltonian H(w) 
depend on a parameter w. Let E{w) be an eigen- 
value of this Hamiltonian. Then the mean value 
of the derivative of the Hamiltonian with respect to 
the parameter can be expressed through the deriva- 
tive of the eigenvalue: {'^ e{w)\-^ H{w)\^ e(w)) = 
-^E{w). Let us now denote the fraction jm as k 
and differentiate the Hamiltonian (2) with respect to 
K. According to the Hellmann-Feynman theorem, we 

get \l^lllldz(d''ldzdz')l,¥{z)il{z')) = dE/dK. 
Now, using {^\z) ^{z')) = (^'''(z - z') ^'(0)), we obtain 
-\ L [{d'^/dz^)gi{z)] = dE/dK, and finally 
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where we have used the known expression for the energy 
(3). 

Note that, as expected, our expressions for the coeffi- 
cients Ci_3 converge, in the limit 7 ^ 00, to the known 
results for the impenetrable bosons (13). This can be eas- 
ily verified using the 7 — > 00 expansion for the function 
e(7) (4). 

Expressions (14), (16), and (15) constitute the second 
principal result obtained in our paper. 

Concluding remarks. Below we present a discussion 
on empirical observation of and applications for the 1/p^ 
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momentum distribution tails, in experiments with har- 
monically trapped atomic gases. 

(a) First of all, we would like to discuss the momentum 
range where the tail should be looked for experimen- 
tally. Relying on the 7 — » oo results [2] and an analysis of 
the predicted-by-Bogoliubov's-theory momentum distri- 
bution (corresponding to 7 ^ 0), we conjecture that in 
the whole range of the interaction strength 7, the high-p 
asymptotics of the momentum distribution corresponds 
to the range of momenta given by p ^ (m/i/?i^)^/^ for 
all 7, where /z is the chemical potential of the system. 

For the case of a harmonically confined gas at 7 ^ 1 this 

leads to 

p » VNpno for 7 > 1 . (17) 

(b) Our zero-temperature results are valid as long as 
the temperature does not exceed the chemical potential, 
ksT <C fj. (for all 7), or, for the Tonks-Giradeau case, 

kBT « Nhoj for 7 > 1 . (18) 

For temperatures comparable to the chemical potential, 
the 1/p^ law should persist within the range (17), but the 
prefactor is not yet known and it is a subject of future 

research. 

(c) Experimentally, the momentum distribution of the 

Tonks-Girardeau gas can be detected either in situ [26] or 
via a ballistic expansion [27]. The latter option requires 
some caution, especially in the Tonks-Giradeau (7 00) 
case. If only the longitudinal confinement is released, the 
interactions will modify the momentum distribution dur- 
ing the expansion, and, as follows from the Bose-Fermi 
mapping [10], the detected momentum distribution will 
correspond to a free Fermi gas. Instead one should release 
both longitudinal and transverse confinements simulta- 
neously. In this case rapid transverse ballistic expansion 
will lead to a quick drop-off of the density, making the 
interactions negligible. 

(d) We believe that the coefficient in front of the high-p 
tail of the momentum distribution (Fig.l) may provide a 
robust experimental identifier of the quantum regime of 
the gas of interest, and, in particular, serve to detect the 
Tonks-Girardeau regime, (i) The high-p tail is not sensi- 
tive to the finite temperature corrections to the correla- 
tion function, which appear predominantly in the low-p 
(long-range) domain, (ii) In experiments with 2D opti- 
cal lattices, where a single cigar-shaped trap is replaced 
by an array of traps [11], the effect of the residual 3D 
mean-field pressure acting during the expansion becomes 
relevant: the high-p part of the momentum distribution is 
far less sensitive to this effect as compared to the low-p 
part, (iii) The theoretical interpretation of the experi- 
mental results is simpler in the high-p case thanks to the 
applicability of the LDA. In the opposite low-p case, the 
LDA leads to entirely wrong predictions [28] . 



Summary. In this paper, we present a short- 
range Taylor expansion (up to the order \z\'^) for the 
zero-temperature correlation function gi{z) of a one- 
dimensional (5-interacting Bose gas (see Eqns. 12, 14, 16, 
and 15). We compute the leading term in the high-p 
asymptotics of the momentum distribution for both free 
(Eqn. 9) and harmonically trapped (Fig.l) atoms. We 
regard the high-p tail of the momentum distribution as 
an efficient tool for identification of the Tonks-Girardeau 
regime in experiments with dilute trapped atomic gases. 
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